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NOTATION 

Power  input 
Thrust 

Tangential  force 
Number  of  blades 
Diameter  (of  rotor  and  shroud) 

Length  of  shroud 
Lengtb'diamoter  ratio  of  shroud 
Hub  radius 

Cylindrical  coordinates 

Nondimensional  cylindrical  coovdiv*tes 

Velocity  of  approach 
Angular  velocity 
Advance  coefficient 
Circulation 

Nondimensional  circulation 

Axi  &l  component  of  induced  velocity 

Radial  component  of  induced  velocity 

Tangential  component  of  induced 
velocity 

Sink  density  per  unit  are>\ 

Cavitation  number 
Drag  coefficient 
Lift  coefficient 
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Drag'll ‘t  ratio 


Thrust  coefficient 


Power  coefficient 


Symbols  for  functions  which  are 
defined  in  the  paper 
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ABSTRACT 

Tho  flow  and  the  forces  which  are  generated  by  a propulsion  system 
consisting  of  rotor,  shroud  and  guido  vanos  are  analyzed.  For  this  purposo, 
tho  components  of  the  system  aro  replaced  by  proper  singularities.  From  the 
component  velocity  fields  and  from  the  characteristic  constants  of  the  singu- 
larities, the  interaction  forces  between  the  components  are  determined  from 
which  the  net  forcos  of  the  unit  follow.  Tho  deduced  expressions  for  thrust 
and  power  input  taken  together  with  the  pressure  increase  at  the  rotor,  which 
arises  from  the  action  of  the  shroud,  and  with  tho  condition  of  cavitation  free 
flow  form  the  basis  for  a method  of  design  of  a propulsion  unit.  To  apply  this 
method,  knowledge  of  both  lift  versus  angle  of  attack  curves  and  of  pressure 
distribution  curves  of  sections  in  cascado  is  necessary. 

Comparison  of  experimental  results  for  the  efficiency  and  for  the  inter- 
action force  between  rotor  and  shroud  aro  in  fair  agreement  with  tho  respective 
analytical  expressions  taking  into  account  the  lack  of  knowledge  relative  to 
the  drag  of  tho  shroud. 


1.  INTRODUCTION 

Recent  interest  by  the  U.S.  Navy  in  the  possible  application  of  shrouded  propellers 
to  various  types  of  naval  vessels  for  the  purpose  of  delaying  cavitation  and  propeller  noise 
has  led  to  a study  of  the  avMlable  theory  of  such  a propulsion  system.  It  was  found  desir- 
able to  consider  further  the  theoretical  aspects  of  this  problem  which  has  resulted  in  the  de- 
velopment of  the  theory  represented  here.  It  is  planned  to  supplement  this  work  with  a pre- 
sentation of  the  theory  as  appliod  to  a specific  design  problem. 

The  forces  on  the  components  of  the  system  and  the  design  data  are  to  bo  determined 
for  a given  net  thrust  of  a propulsion  system  consisting  of  rotor,  shroud,  and  guide  vanos 
(the  latter  two  components  being  stationary)  and  for  given  quantities  of  speed  of  advanco, 
number  of  revolutions,  rotor  diameter,  and  pressure  at  the  rotor  (duo  to  the  action  of  the 
shroud).  For  this  problem,  it  is  necessary  to  ascertain  the  mutual  interaction  betwoen  the 
components  of  the  system  which  follows  when  those  components,  with  respect  to  thair  effect 
on  the  flow,  are  replaced  by  proper  singularities.  In  order  t j deduco  the  component  forces  on 
the  basis  or  such  a theory,  it  is  necessary  to  know  the  component  velocity  fields.  This  re- 
quires certain  approximations  which  will  be  mentioned  first  i'.  connection  with  the  singulari- 
ties. 

The  shroud  is  considered  to  he  an  annular-shaped  thin  hydrofoil;  all  tho  sections, 
determined  by  a meridian  cut,  are  supposed  to  be  of  equal  shape.  This  axisymmetrical  hydro- 
foil is  replaced  by  a tow  of  adjacent  circular  vortices.  Furthermore,  it  is  assumed  that  the 
shroud  is  a circular  cylindor  with  a diamoter  equal  to  that  of  the  rotor.  This  approximation 
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holds  if  both  the  ordinates  of  the  camber  lim?  of  the  hydrofoil  section  as  measured  from  the 
cylinder  and  the  square  of  the  slope  of  the  ordinates  are  small.  Then,  the  circular  vortices 
can  be  arranged  on  the  surface  of  the  cylinder,  i e.,  their  diameter  is  constant  and  equals 
that  of  the  rotor. 

The  gap  between  rotor  and  shroud  is  assumed  to  be  very  small  and  the  length  of  the 
shroud  to  be  sufficiently  great  so  that  a radial  flow  around  the  tips  of  the  rotor  blados  and 
the  equalizing  of  pressuro  which  arises  from  such  radial  flow  are  prevented.  Then  a circu- 
lation which  is  independent  of  the  radius  is  po&siblo  at  the  rotor  for  which  free  vortices  ore 
not  present.  When  the  blades  are  replaced  by  lifting  lines,  the  vortex  system  of  the  rotor 
then  consists  merely  of  the  lifting  lines  and  of  a hub  vortex  of  the  combined  circulation  of 
that  of  the  lifting  lines.  Further,  discontinuities  of  the  velocity  components  occur  within 
the  boundary  of  the  rear  part  of  the  jet.  The  introduction  of  helical  vortex  lines  within  the 
boundary  would  be  necossary  to  account  for  theso  discontinuities. 

With  respect  to  the  velocity  Held  of  the  rotor,  the  influence  of  a finite  number  of  blades 
is  neglected;  by  this  simplification,  the  field  of  the  absolute  velocity  becomes  independent 
of  time.  It  is  shown  later  that  this  assumption  may  be  justified  with  the  number  of  blades  in 
practical  application.  Then  the  discontiuuity  cf  the  velocity  within  the  be  -ndary  of  the  roar 
jet  is  equivalent  to  a vortex  sheet  consisting  of  helical  vortex  lines.  It  is  ■'**>,1  known,  that 
such  a sheet  can  be  resolved  into  two  sheets,  one  consisting  of  straight  vorte>  ices  which 
are  parallel  to  the  axis  and  the  other  one  of  ring  vortices  which  are  perpendicular  to  the  axis. 

The  first  sheet  produces  only  tangential  velocity  components,  as  follows  from  she  law 
of  Biot-Savarti  Within  the  inside  fluid,  i.e.,  within  the  slipstream,  the  velocity  from  tins 
sheet  equals  zero  (as  follows  from  Stokes’  law).  The  tangential  velocity  components  pu 
ducsd  from  this  sheet  on  the  shroud  (which  is  within  the  outside  fluid)  are  without  interest 
for  the  interaction  between  rotor  and  shroud.  Hes.ce,  the  sheet  consisting  of  straight  vortex 
lines  need  not  be  considered  in  the  following. 

The  sheet  of  ring  vortices  can  be  replaced  by  singularities  on  the  propeller  disk.  It 
is  well  known  that  the  velocity  potential  of  a single  dosed  vortex  equals  that  of  a uniform 
distribution  of  dipoles  over  any  surface  bounded  by  the  vortex,  the  axes  of  the  dipoles  being 
perpendicular  to  the  surface  (e.g.t  see  Lamb,  Hydrodynamics,  Art.  150,  161,  and  102).  The 
singularities  on  the  disk  which  correspond  to  a row  of  dosoly  adjacent  ring  vortices  behind 
the  disk  follow  when  integrating  the  potential  of  a single  ring  vortex  from  g = 0 (at  tho  rotor 
disk)  to  infinity  downstream,  as  given  in  Lamb,  Art.  161  (13).  The  result  for  the  inflow  is 
that  the  potential  of  this  tow  of  ring  vortices  equals  that  of  a uniform  sink  distribution  over 
the  disk.  The  same  result  has  been  obtained  by  Dickroan  in  a different  way,  viz.,  from 
Euler's  differential  equations. 1 
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In  the  following  discussion,  the  rotor  inflow  is  br'sod  on  the  flow  of  a sink  disk  by 
which  numerical  calculations  bocomo  simpler  thnn  from  the  velocity  field  of  the  cylindrical 
vortex  sheet  at  the  boundary  of  the  roar  jot,  Tho  forco  on  the  rotor, however,  is  not  calcu- 
lated from  the  forces  which  acton  tho  sinks  but  from  thoso  on  tho  lifting  vortices,  permitting 
the  influence  of  a finite  number  of  blades  on  tho  forco  to  bo  approximated. 

Analogously,  the  guide  vanos  oro  replaced  by  lifting  lines  in  order  to  determine  the 
forces.  The  total  circulation  of  those  lifting  linos  is  taken  oqual  to  tnat  of  tho  rotor  in 
which  case,  with  a finite  number  of  blades,  tho  average  of  the  tangential  velocity  component 
within  tho  slipstream  equals  zero  at  tho  design  advance  coefficient.  With  respect  to  the 
influence  of  the  guide  vanes  on  the  axial  and  radial  inflow,  which  depends  on  the  thrust 
loading,  rotor  and  vanes  are  considered  as  a unit  which  causes  tho  sum  of  the  component 
thrusts.  Correspondingly,  tho  strength  of  tho  sinks  at  the  rotor  disk  is  determined  from  the 
sum  of  tho  thrusts  of  rotor  and  vanos.  This  amounts  to  a nogloct  of  both  the  finite  distance 
between  rotor  and  vanes  and  of  the  variation  of  tho  circulation  with  time  which  arises  from 
the  interference  betwoon  rotor  and  vanos  with  a finite  number  of  blades. 

The  vanes  are  considered  to  be  situated  behind  tho  rotor.  In  this  position,  the  vanes 
cause  a thrust  which  is  in  tho  same  direction  as  tho  thrust  generated  by  the  rotor.  On  the 
other  hand,  due  to  the  pressure  increase  which  arises  from  tho  positivo  thrust,  the  losses 
at  the  vanes  are  greater  than  when  located  in  front  of  tho  rotor.  In  this  latter  arrangement, 
a negative  thrust  arises  at  tho  vanes  which  leads  to  a pressuro  drop  and,  therefore,  to  small 
losses.  At  the  same  time,  the  necessary  lift  coefficients  at  the  rotor  are  smaller  in  the 
latter  arrangement  bocause  of  a greater  relative  velocity;  this  advantage  with  respect  to  the 
onset  of  cavitation  may  be  offset  by  the  smaller  cavitation  numbers  of  tho  sections  of  the 
rotor. 

The  termination  of  the  rotor  shaft  requires  singularities  which  are  chosen  so  that 
the  meridian  cut  of  tho  shaft  becomes  a streamline  of  the  combined  flow.  If  this  condition 
is  not  satisfied,  tho  continuity  equation  is  violated.  Approximately,  however,  a single  sink 
on  tho  axis  may  be  sufficient  to  account  for  the  effect  of  tho  shaft. 

When  the  components  of  the  system  are  replaced  by  the  afore-mentioned  singularities 
within  a uniform  velocity  field,  which  is  identical  with  the  speed  of  advance,  the  force  act- 
ing  on  any  one  singularity  follows  from  a general  rule  (Lagally).  For  the  problem  under 
consideration,  the  main  consequence  of  this  general  rule  is  that  the  interference  velocity- 
on  each  of  the  singularities  arising  from  tho  other  singularities  present  does  not  lead  to  a 
resultant  force  but  to  interaction  forces  between  the  singularities.  From  this,  the  reaction 
between  shroud  and  rotor,  e.g.,  is  determined  either  from  the  singularities  of  tho  shroud  and 
the  velocity  induced  at  the  shroud  from  tho  rotor  or  from  the  singularities  of  the  rotor  and 
the  velocity  induced  at  tho  rotor  from  tho  shroud.  Tho  force  follows  from  the  product  of  the 
characteristic  constant  of  the  singularity  (e  or  I’)  times  the  velocity.  In  the  case  of  a sink 
singularity,  the  direction  of  the  force  equals  that  of  the  volocity;  in  the  case  of  a vortex, 
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the  forco  is  normal  to  the  velocity  and  coincides  with  the  vector  product  (»  x curl  v).  This 
latter  special  case  of  the  general  rule  is  the  well-known  law  of  Kutta-Joukowsky. 

l THE  VELOCITY  FIELDS 

2.1  THE  VELOCITY  FIELD  OF  THE  SHROUD 

The  velocity  field  of  the  shroud  is  determined  from  its  chordwise  circulation  distribu- 
tion; two  esses  can  be  considered  in  this  respect,  viz.,  the  shape  of  the  ahieud  is  given  and 
the  circulation  distribution  is  to  bo  determined  or  vice  versa.  In  both  of  these  cases,  the 
unknown  quantity  follows  from  the  boundary  condition  of  the  flow,  viz.,  that  the  section  of 
the  shroud  is  a streamline  of  the  relative  flow  at  the  shroud.  This  Row  results  from  the 
undisturbed  flow,  from  the  induction  of  rotor,  guide  vanes,  and  shaft  sink,  and  from  the  self- 
induction  of  the  shroud. 

The  problem  of  determining  the  circulation  distribution,  neglecting  the  thickness  of 
the  shroud  section  end  considering  the  camber  line  given,  leads  to  a complicated  integral 
c^uaucn  Cc*  which  Oickmann  gives  an  approximate  ? vlution.2  However,  teas  solution  be- 
cc'r.o a very  laborious  when  h « 1/2R  £ 1,  h being  the  ratio  of  shroud  leng  h to  its  diameter. 
This  case  cannot  be  excluded  from  consideration.  Therefore,  it  ir,  desirable  to  avoid  begin- 
ning with  a given  camber  line  and  to  begin  instead  with  a »ivei  circulation  distribution 
which  necessitates  the  determination  of  the  camber  line  altv  anas  such  that  the  circulation 
distribution  is  realized,  i.o.,  that  the  afore-mentioned  boundary  condition  of  the  resultant 
flow  is  satisfied.  Proceeding  in  this  way  avoids  the  integral  equation,  and  the  camber  line 
follows  from  a first-order  differential  equation. 

The  circulation  function  over  the  chord  length  of  the  shroud  is  chosen  to  be  a half 
ellipse  in  order  to  avoid  peaks  of  negative  pressure.  This  circt«:  tion  function  requires 
that  the  forward  stagnation  point  be  situated  at  the  leading  edge  of  the  camber  line  (shock* 
flee  condition).  It  is  well  known  that  this  condition  is  satisfied  in  the  case  of  two-dimensional 
flow  at  a flat  wing  with  a parabolic  camber  line  at  zero  angle  of  attack.  The  camber  line  and 
the  geometric  angle  of  attack  of  a ring-shaped  wing  which  satisfy  the  sbock-froe  condition 
are  not  yet  known.  Both  of  those  unknowns  can  be  determined  after  the  component  velocity 
fields  have  been  established. 

For  the  velocity  field  of  the  shroud,  a cylindrical  sheet  of  adjacent  circular  vortices 
is  considered  corresponding  to  the  afore-mentionBd  assumptions.  Tho  stream  function  of  a 
single  circular  vortex  of  radius  ft,  width  da',  and  circulation  (yd a')  is,  in  cylindrical  coor- 
dinates, r,  a,  and  <p,  represented  by  (e.g.,  seo  Lamb,  Art.  161) 

Tr-(yda')(rRV  k(ZD-K)  with  D = 

— tit  I.* 
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K and  E are  the  complete  elliptic  integrals  of  the  first  and  second  kind,  respectively, 
of  modulus  k which  is  determined  by 

*rR 

=(a-a‘)*  + (r  + K)2 

where  a ' fixes  the  location  of.tho  vortex  ring. 

From  this  expression  for  the  stream  function,  tho  velocity  components  of  a single 
circular  vortex  ring  are  deduced  when  applying  the  following  relations  for  the  derivatives 
of  the  complete  elliptic  integrals  with  rospect  to  the  modulus 

d(kD)  _ d^E  K-D 
dk  ““  dk2  ” 1-k2 


d(fcK)  E 
dk  l-k'- 


One  obtains 


dwa  - t |2rD+ (R-r) 

. 2*21  f g _ » Ti\ 

” t(o— a')2  + (r  +R  )Z]F  i'1”***  ' 


dwt  “0 


Because  of  the  axial  symmetry  of  the  problem,  both  the  axial  and  radial  components 
ate  independent  of  the  angular  coordinate,  and,  further,  tho  tangential  component  equals 
zero. 

For  a vortex  cylinder  of  any  circulation  distribution,  the  velocity  components  follow 
from  the  elementary  components  of  a vortex  ring  by  integration.  The  result  written  nondimen* 
sionally  for  an  elliptic  distribution  is; 
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x 

where 

z 


are  nondimensiontl  cylindrical  coordinates, 


l 

k* 


is  the  length  of  the  shroud, 

4 % 


18 


[*(*-«')]  + (l+*)‘ 
is  1 - A2, 


h - 


2 R 


* " 1/2 

rs 

Gs“  ( 2Rv .) 


is  the  ratio  of  length  to  diameter  of  the  shroud, 

is  the  variable  of  integration  over  the  chord  length  of  the  shroud,  and 

is  the  nondimensional  total  circulation  of  the  shroud.  This  quantity  is 
positive  if  the  force  which  is  generated  by  a positive  vQ  (in  the  direction 
of  positive  a)  is  positive  (in  the  direction  of  positive  r). 


The  complete  elliptic  integrals  within  the  velocity  components  have  been  expressed 
by  integrals  over  a product  of  the  modified  first  order  Bessel  functions  in  a paper  by  Stewart.3 
For  numerical  purposes,  however,  the  velocity  components  as  written  in  Equations  [1]  and 
[2]  are  found  to  be  more  convenient. 

Relative  to  the  problem  under  consideration,  the  velocity  components  of  interest  are 
those  within  the  plane  normal  to  the  axis  through  the  halfway  point  of  the  shroud  length, 
a - 0,  (in  which  plane  rotor  and  vanes  are  supposed  to  be  situated)  and,  further,  at  the  vor- 
tex sheet,  * - 1. 

For  the  plane  a - 0,  there  is  obtained  from  Equations  [1]  and  [2]: 


-£/>’ 2 l/Tr?T  + 
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The  radial  component  in  this  plane  becomes  zero  since  the  integrand  is  an  odd  func- 
tion of  s'. 

At  the  axis,  x ■»  0,  the  integrand  for  wa  can  be  oxpressod  by  known  functions. 

One  obtains 


o 
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Vl-z'd* 

[(I^“ 

yT+A* 
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A*)  E(k*)  - 


At  stations  x which  are  different  from  zero,  the  integral  for  wa  can  only  be  evaluated 
by  numerical  calculations.  The  results  of  these  calculations  over  a range  of  A are  represent- 
ed on  Figure  2.  From  this  figure  it  follows  that  the  velocity  distribution  within  the  middle 
plane  becomes  uniform  for  shrouds  of  a length  greater  than  about  twice  the  diameter.  For 
shorter  shrouds,  the  axial  velocity  becomes  increasingly  less  uniform  the  shorter  the  shroud 
is,  tho  volocity  increasing  from  the  axis  towards  the  wall  of  the  shroud. 

It  should  be  mentioned  that  different  notations  for  the  same  velocity  component  are 
used  in  this  report.  For  instance,  the  afore-deduced  velocity  component  0 is  the 


Figure  1 Figure  2 * Velocity  Field  of  the  Shroud  in 

the  Plane  a - 0 

Middle  plane  normal  to  the  exla. 
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component  induced  at  the  plane  3 - 0, 4.9., 1 at  the  disk  fn-m  the  shroud.  In  the  application 
of  this  work,  a notation  with  two  indices  is  found  more  convenient.  The  first  index  fixing 
the  plane  of  reference  and  the  second  one  the  component  of  the  system  by  which  the  respec- 
tive velocity  component  is  generated.  Correspondingly,  the  velocity  component  («>„),„(,  is 
denoted  on  the  diagrams  and  later  on  in  this  report  by  (toa)p,s  Analogously,  (tca)s  s denotes 
the  self-induction  of  the  shroud, 0 the  axial  component  induced  at  the  shroud  from  the 
disk,  and  so  on. 

The  self-induction  of  the  shroud  follows  from  Equations  II]  and  [2]  with  * « 1 


= — l/ 1 — z'2  kaD(k)  d( z‘—  z) 

a*!  y - - 2o,»} " 


k2  m i 

[*u-«*,r+4 


Both  of  these  integrals  become  improper  integrals  when  3' approaches  3.  In  this  case, 
k approaches  1,  A' approaches  zero,  and  K approaches  infinity  as  nat  log  [1/(3  - s')].  With 
respect  to  the  integral  for  to,,,  the  transformation  («'-  3)  - l*  is  introduced.  Then  the  product 
(K  t*~ l)  approaches  zero  when  a'  approaches  3 provided  that  n > 1. 

Correspondingly,  n = 3 is  chosen  for  the  numerical  evaluation  of  the  integral  in  order 
to  have  real  values  of  t for  negative  quantities  (s'  - 3);  it  follows  that 

\—  7f'  ” " ^ J », KT— z'*A*D(fc)  t2dt 

vu  “si  = 1 " /iT* 

The  results  of  the  numerical  evaluation  are  represented  in  Figuro  8.  The  curves  are 
symmetrical  about  a - 0. 

Rolativo  to  to,,  the  integral 

//  - f+l  (s'-  zWl-z'2  k3  dz' 

J~l  (kf>* 

is  considered  first;  the  integrand  oecomes  indefinite  of  form  0/0  when  s'  approaches  3. 
Replacing  { k ')2  by  (l  - k2)  gives 


H 


4 

A7 


f ks  E 


♦ 1 


dz  '+ 


r- 1 


z — z 


J-J 


(z'-z)^l  -z'2  ka  E dz ■ 
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Figure  3 • Self-Induction  of  the  Shroud 

Axial  Compcnant 


The  First  of  thece  two  integrals  is  an  improper  integral.  To  determine  its  principal 
value,  we  introduce  a'-  cos  <p\  da'~  - sin  <p dip',  a = cos  and  resolve  the  numerator 
(yjl  - a‘2E k3 aintfi)  within  the  interval,  U>'»  n to  V>'-  0,  which  corresponds  to  tho  interval 
a»-lto2«  + l,  into  a Fourier-series 

. > 
y 1 — *'*  k*E  ain  tp  =>  Jfc*E  ain*^'  = y^gHcos  rup 

*1  3 0 

Then  it  follows  that 


r 


con  ny 


Jo  cos  ip  — cos  <p 
With  this  expression,  there  is  obtained 

’ Wr  1 


«-Z 


_ _ am  ny 

3a  " » ^ 

TTF-o  «««♦» 


f-Xa,  si; 


am  ne 


( ? ^ ^ * w *0  ‘t* 2 + ‘!>b  • “1  *} 


Within  the  integral,  tho  product  A'(a'-  a)  approaches  aero  when  a' approaches  a;  thus 
the  integral  can  be  evaluated  numerically  without  difficulty.  Relative  to  the  first  term  within 
the  brackets,  this  term  becomes  indefinite  of  form  0/0  at  the  end  points  of  the  interval 
a o - 1 and  a - + S.  From  t’lloapitai’s  rule  it  follows  that 


CONFIDENTIAL 


CONFIDENTIAL 


10 


for  z =■  - 1: 


for  z = + 1: 


The  numerical  results  as  obtained  by  means  of  Fourier-expansions  consisting  of  six 
terms  an  aro  represented  on  Figure  4.  Since  the  radial  component  is  an  odd  function  of  a, 
it  is  represented  only  for  an  interval  of  3 from  - 1 to  0. 

It  should  be  mentioned  that  the  main  difference  between  a ring-shapod  and  a flat  wing 
lies  in  tho  axial  component  of  the  self-inducod  flow.  This  component  is  zoro  in  the  case  of 
a flat  wing  and  arises  from  the  lateral  parts  of  the  ring  for  a ring-shapod  wing. 

The  velocity  which  is  induced  at  the  shaft  sink  from  the  3hroud  is  determined  approxi- 
mately by  assuming  the  total  circulation  of  tho  shroud  I'  concentrated  in  a *inglo  vortex  at 
3 = 0.  From  the  stream  function  of  a singlo  circular  vortox,  one  obtains  at  tho  shaft  sink 
with  coordinates  xsh  = 0 and  zsh°  ath! 


(!2s!\ 

'"o<V,a,s 


f la] 


2.2  ThE  VELOCITY  FIELD  OF  ROTOR 
AND  GUIDE  VANES 

As  mentioned  previously,  the  effect 
of  rotor  and  guide  vanes  on  the  inflow  is 
deduced  from  that  of  an  axisymmetrical  sink 
distribution  over  tho  rotor  disk.  For  a disk 
of  small  thickness  8,  the  volume  element 
amounts  to  dF  »=  r'(dr')  (dtp)  8. 

Lot  g be  the  strength  of  the  sinks 
per  unit  volume;  then  the  strength  of  the 
volume  element  becomos  de  = r'(dr')  (dtp) 
(q8).  basing  to  the  limit  8 = 0 in  such  a 
way  that  the  “surface  density”  e ^ = q8 
remains  constant,  the  elementary  strongth 
of  the  surface  element  d A = r'(dr')  (dtp) 
becomes  (eD  d A)  and  induces,  at  a point 
with  cylindrical  coordinates  (a,  r,  0),  the 
velocity  element 


Figure  4 - Self-Induction  of  the  Shroud 


Redial  Component 
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(d  w) 


cp  d A * t 


4np * 4np' 


'■}  r dr  'dtfi 


where  p is  the  distance  betweon  the  area  olement  of  the  disk  and  the  point  of  reference.  If 
the  area  element  has  coordinates  (0,  r\ <p)>  it  follows  that 

2 v 2 

p*  = a2  + ( r'ain  <fi)  + ( >’  - x'coe  <p) 

For  the  components  of  (d2tc),  there  is  obtained: 


(d*w)a  = (d*w)  - 
P 


(rf2w),=  (d2u>)  cos  (r‘,p)  = (d*w)  (--  £ cos<p) 


Introduction  of  the  expressions  for  (d2tc)  and  p into  those  relations  and  integration 
over  <p  yield  the  following  result  for  the  velocity  Held  of  an  annular  sink  ring,  eD  being  in* 
dependent  of  <p 


2 Uz)(x)2 


Elk) 


dW*~^dXTh*?+(x--x?  ^+(xVx)*  2irx 
dwr  = - P-.  dx 


e°—  dxw; 


2**  Vuizf+ix'+x)2 


{ 


K(k)  - E(k) 


[- 


2x'(x'-x) 


]}  - ix'w' 


[31 


[41 


{hzf  + (x’-  x)2 
The  modulus  k of  the  complete  elliptic  integrals  has  the  value 

k _ 2V/xxr 

^(Ax)*+  (x’+x)2 

The  tangential  component  is  zero  because  of  axial  symmetry. 

From  these  expressions,  the  velocity  components  of  a sink  disk  follow  when  integrat- 
ing over  x'.  Since  eD  may  depend  on  x',  there  is  obtained 


w* 


2*J, 


W. 


2rrJ: 


*h 


f-dz 

[3al 

Hdz 

X 

[4al 

where  the  nondimensional  radius  of  the  hub  th/R  is  denoted  by  xA, 

Putting  z » 1 in  the  expressions  for  it*  and  it*  gives  the  velocity  components  which 
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Figure  5 

are  induced  at  the  shroud  from  the  disk.  For  this  case  and  for  A » 1,0,  the  one  factor  of  the 
integrands,  viz.,  {v*/z‘)S  D and  (v>*/z'}S  D is  represented  on  Figures  5 and  6,  respectively. 
In  these  representations,  most  of  the  numerical  values  ot  the  functions  tc*  and  tc*  have  been 
taken  from  a paper  by  Kuechemann.4  With  respect  to  the  influence  of  A,  the  determination  of 
the  curves  for  one  value  of  A is  sufficient  since  te*  and  tr*,  from  Equations  [3]  and  [4] , re- 
spectively, depend  merely  on  the  product  (A  a).  That  is,  for  different  valuos  of  A,  the  curves 
for  u * and  te*  shift  so  that  the  same  ordinate  belongs  to  a different  abscissa,  which  follows 
from  (A  a)  « constant. 

In  order  to  determine  («0)S(£,  and  («,)s,o*  the  variation  of  eD  over  the  radius  must  be 
known.  If  eD  is  independent  of  aroi  can  bo  replaced  by  a suitable  average  eD,  the  expres- 
sions become 
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'w*Ks,i,jr 


[3b! 


, , 2n 

et> 
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which  can  be  computed  numerically  from  the  quantities  on  Figures  5 and  6;  the  rosult  is 
given  in  Figures  7 and  8,  respectively,  for  a hub  radius  xk  = 0.4  and  for  a range  of  h.  With 
respect  to  this  latter  variable,  the  curves  are  related  in  the  some  manner  as  the  integrands. 

Relative  to  the  aolf-iaductio"  of  the  sink  disk,  the  axial  component  of  the  induced 
velocity  is  of  interest  because  of  its  relation  to  tho  thrust.  From  Equation  [3J,  the  compo- 
nent would  follow  at  each  station  x from  : 

Instead  of  determining  this  limit  of  tho  integral,  thoro  is  a moro  convenient  wr.y  of 
ascertaining  „0  by  representing  wa  by  2 series  of  Logoudro  polynomials' 


coNHBftrruu 
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Figure  7 
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For  ft  sink  disk  of  radius  It  with  c{)  assumod  to  bo  indopondont  of  tho  radius,  the 
volocity  noton  tin!  at  any  point  r<  It  is  reprosontod  in  spherical  coordinates  r,  «>  ,<p  by 

tho  sorios 


4> 


<r,t» 


1 r ‘ 1 r* 

2 /PP*~  gR*P*  + 


t’n  « V n (cos^)  is  tho  Logondro  polynomial  of  n,h  dogroo.  This  potontial  satisfies  Laplaco's 
equation,  tho  boundary  conditions  boing  tho  values  of  tho  potontial  at  infinity  and  at  tho 
axis;  tho  latter  is  known  from  a simple  integration  ovor  tho  radius  of  tho  disk.  By  comparison, 
it  is  seen  that  thoso  boundary  conditions  arc  identical  with  thoso  of  tho  gravitation  potontial 
of  tho  disk;  tho  dotails  of  this  deduction  aro  available  in  standard  toxtbooks. 

From  tho  aforo- written  sorios  for  <*>  (which  is  indopondont  of  <fi  and,  thoroforo,  axi- 
symmotrical),  tho  axial  volocity  component  at  points  arbitrarily  noar  tho  disk 
(t>  =•  «,«•*>  0)  becomes 

. ! . . = \ for  r < * 15] 

t 

• * o 


For  points  of  roforonco  for  which  r > It,  tho  potential  is  roprosontod  by 


a _ tu  Rt  i a _ i ,, . i *’ \ 

0(r,i»  2^*2  r 8ri/*  + I«7r,« ) 


from  which  thcro  is  obtoinod  for  points  arbitrarily  noar  tho  plane  of  tho  disk 

(».)/.,«  a *0  for  r > It 


t61 


If  Cjj  is  indopondont  of  tho  radius,  tho  axial  volocity  componont  is  constant  ovor  the 
disk  and  is  zero  outside  of  tho  disk  for  points  within  its  piano. 

Now  compare  the  volocity  field  of  a sink  disk  with  that  of  a rotor.  At  tho  disk,  tho 
velocity  jumps  from  -(c^)/2  to  »(e0)/2;  this  follows  from  tho  aforo-writton  volocity  poten- 
tials. Otherwise,  the  volocity  fiold  of  tho  sink  disk  is  continuous.  Jn  particular,  with  re- 
poet  to  tho  axial  componont,  it  follows  from  Equation  (3]  that  positive  values  of  z (which 
aro  behind  the  disk)  yield  a negative  quantity  (tc0),  eD  being  a sink  donsity.  As  compared 
to  the  absolute  flow  of  a propeller,  this  holds  for  the  volocity  field  behind  tho  propeller 
outside  of  tho  slipstroam  but  not  within  tho  slipstream.  To  approximate  tho  slipstream  bo- 
hind  the  rotor  with  respect  to  tho  axial  component,  it  is  necessary  to  suporimpose  on  tho 
sink  flow  behind  tho  disk  a volocity  field  of  magnitude  |e/;|  which  is  in  tho  direction  of 
positive  2.  Then  tho  axial  volocity  componont  right  behind  tho  disk  bocomes 
(|«0|  - jc0j/2)  = |e^|/2  which  equals  tho  inflow  into  tho  disk.  Infinitely  far  behind  tho 
disk,  the  effect  from  the  sinks  is  ze.o  and  tho  volocity  equals  \?p\.  Hence,  tho  axial 
volocity  of  the  combined  flow  behaves  like  tho  axial  velocity  componont  of  a propollor,  viz., 
it  is  continuous  at  the  propollor  and  infinitely  aft  it  increases  to  twice  the  valuo  at  tho 
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propeller.  Further,  considering  the  axial  component  ii  the  vicinity  of  a cylinder  of  radius 
/?,  the  axial  component  of  the  absolute  flow  is  continuous  when  traveling  through  the  cylin- 
der  in  front  of  the  disk  but  discontinuous  behind  it.  This  complies  with  the  properties  of 
the  axial  velocity  component  of  a propeller. 

Behind  the  disk,  the  axial  component  of  the  absolute  propeller  flow  jumps  from 
[Ifot  + wa(  + s)] inside  to  ti>0(+  3)  outside  of  the  raco,  and  at  the  boundary  it  equals 
[<*[)/ 2 *■  te0(+3)].  Since  tcj.+  a)  « - ua(-  a)  the  axial  component  of  tho  interference  velocity 
at  the  shroud  from  the  rotor  for  positive  a expressed  in  terms  of  the  sink  flow  becomes 


wa(~  z )*,/< 


[3c] 


tc0(-  a)s  0 following  from  Figure  7. 

in  addition  to  the  axial  and  radial  velocity  components,  there  is  also  a tangential 
velocity  Field  from  rotor  and  guide  vanes  which  is  caused  by  the  hub  vortices.  When  rotor 
and  guide  vanes  are  considered  separately  and  the  blades  are  replaced  by  lifting  lines,  the 
following  conclusions  are  known  from  Stokes'  Law  with  respect  to  the  tangential  velocity 
field  of  each  of  the  components:  The  field  is  zero  in  front  of  the  lifting  lines  and  is 


wt  m 


nr 

4rn 


[7] 


at  a lifting  line,  where  r is  the  circulation  at  one  of  the  lifting  linos  and  n is  the  number  of 
blades  on  rotor  or  number  of  vanes.  In  the  flow  behind  the  respective  component,  the  tan- 


gential velocity  equals  twice  this  value. 

The  velocity  which  is  induced  at  the 
shaft  sink  from  the  sink  disk  follows  from 
Equations  [3a]  or  [3b]  when  introducing 

z . » « ii*  _ 

**  lf2  1/2  2/i 

and  when  adding  the  (negative)  quantity  so 
obtained  to  the  uniform  velocity  Hold  |eD|. 
(Relative  to  the  position  of  the  shaft  sink  on 
the  axis,  see  Section  2.3.) 

When  a sink  density  at  the  disk  which 
is  independent  of  the  radius  cD  is  assumed, 
the  nocdimensional  axial  component  at  the 
shaft  sink  becomes 


Figure  9 
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This  function  F is  represented  on  Figure  9 for  a ra.igo  of  tho  ratio  ^jA/(i/2).  i.o,, 
the  ratio  of  the  shaft  length  measuiod  from  2 - 0 to  Uto  length  of  half  tho  shroud,  and  for 
several  values  of  A,  being  chosen  0.4. 

2.3  THE  VELOCITY  FIELD  OF  THE  SHAFT  SINK 

The  termination  of  the  shaft  requires  a sink-source  distribution  on  the  axis.  Whon 
this  distribution  is  approximately  replaced  by  a single  sink,  it  follows  frr  m well-known 
relations,  zh  « K tiK/H  being  the  nondimensional  radius  of  the  shaft,  that 


for  the  strength  of  the  sink 

= R'^x* 

[8] 

for  the  position  of  tho  sink 

, ..  2*  m La  _ A ft 

<*  i/2  1/ 2 2 A 

[9] 

These  relations  hold  whon  the  velocity  vfli  is  uniform.  With  the  rotor  and  shroud 
working,  tho  difficulty  arises  that  the  velocity  field  within  which  the  shaft  is  situated  is  no 
longer  uniform  but  depends  on  both  the  radial  and  axial  coordinates.  As  an  approximation, 
the  sum  of  the  velocity  of  approach  and  of  tho  axial  components  which  are  induced  at  the 
place  of  the  shaft  sink  from  shroud  and  rotor  will  bo  introduced  for 


!V»  4.  j _ G.-  + JL  If*  F 

vo  ' (l  + *,2AA2)f  2w  k ' 


[10] 


the  function  F t being  represented  on  Figuro  9. 

This  approximation  does  not  satisfy  the  boundary  condition  that  the  normal  component 
of  the  flow  relative  to  the  shaft  is  zoro.  Tho  errors  which  arise  should  bo  checked  after- 
wards by  calculating  the  normal  component  of  the  combined  flow  on  tho  given  shaft  contour. 

With  these  relations  for  the  shaft  sink,  one  obtains  for  the  velocity  components  which 
aro  induced  at  the  shroud  from  tho  shaft  sink 
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The  foregoing  considerations  and  numerical  evaluations  give  a sufficient  knowledge 
of  the  velocity  fiold  of  the  shrouded  propeller.  This  velocity  field  forms  tho  basis  for  ascer- 
taining the  force  components. 
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3.  THE  FORCES 

As  mentioned  on  page  3,  the  interaction  forces  result  from  the  mutual  interference 
velocities  botween  the  singularities.  The  interactions  with  the  shaft  sink  will  be  omitted 
whon  determining  these  forces  since  these  interactions  are  very  small.  Furthermore,  these 
quantities  cannot  bo  determined  very  accurately  because  of  the  approximations  for  the  strength 
of  the  sink.  (The  velocity  Held  of  the  shaft  sink  will,  however,  bo  taken  inu>  account  for  the 
camber  line  of  the  shroud  since  in  this  case  it  is  essential  in  order  to  satisfy  the  continuity 
equation.) 

Disregarding  the  interactions  with  the  shaft  sink,  the  net  thrust  of  the  system  will  be 
obtained  when  the  axial  force  on  the  shroud  is  subtracted  from  the  thrust  of  rotor  and  guide 
vanes. 


3.1  THE  FORCES  ACTING  OH  THE  SHROUD 

The  axial  force  on  the  shroud  vortices  Ts  arises  from  the  radial  velocity  component 
which  is  induced  at  the  shroud  from  other  singularities  present. 

In  this  case,  the  singularities  present  are  the  oinks  at  the  disk  (which  replace  rotor 
and  guide  vanes)  and  the  hub  vortices.  Since  the  direction  of  the  velocity  induced  from  the 
latter  coincides  with  the  direction  of  the  vortex  lines  of  the  shroud,  the  effect  on  the  force 
is  zero.  Then  with  the  radial  component  from  the  sink  disk  (tf,)5>/) , the  axial  force  becomes 

Ts,u  “ 2pH2nhj  riz'AwJ^ds  (12a) 


With  a positive  y(z)  (which  represents  the  circulation  per  unit  length  over  the  shroud] 
and  with  an  inward  radial  velocity  component  from  the  sink  disk,  the  axial  force  on  tho  shroud 
is  directed  backwards,  i.e.,  this  force  represents  a resistance. 

Introducing  for  y(a)  the  elliptical  distribution  which  is  assumed  in  this  paper,  viz., 


y(r)  - 


one  obtains  for  the  nondimensional  coefficient  of  the  axial  force  on  the  shroud  which  is  in* 
duced  from  the  disk 
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An  alternative  expression  for  this  force  can  be  deduced  when  considering  the  reaction 
of  the  sink  disk  to  the  interference  velocity  from  the  shroud  vortices.  From  the  remarks  on 
page  3 , this  reaction  force  is  equal  to  but  opposite  to  that  which  is  caused  by  the  interference 
velocity  from  the  sink  disk  at  the  shroud  vortices  and  which  is  expressed  by  Equation  [12a]. 
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For  this  reaction  force  (denoted  by  TD  s)  in  terms  of  quantities  of  the  disk,  th  > following  is 
obtained 


[13a] 


[13b] 


K«fW -K*>* 


[14] 


Both  of  the  expressions  [12b]  and  [13b]  depend  on  the  radial  distribution  of  the  sinks 
over  the  disk,  Equation  [13b]  in  a direct  way  and  Equation  [12b]  by  (tc,)s  D which,  from 
Equation  [4a],  is  rotated  to  <0(z).  In  order  to  evaluate  Equations  [12b]  and  [13b],  the  distri- 
bution eD(x)  must  be  known.  In  the  case  that  eD  is  independent  of  x or  can  be  replaced  by 
an  average  eD,  the  force  coefficients  can  be  evaluated  once  and  for  all.  In  this  case  if  fol- 
lows that 


(cr)s 


v0 


•v  Gs 


l 


[12c] 


and  that 


[13c] 


The  factors  in  brackets  are  known  from  Figure  8 and  Figure  2,  respectively.  Tho  representa- 
tion of  these  expressions  on  Figure  10  shows  that  Equation  [14]  holds  within  the  accuracy 
of  numerical  integration  which  is  limited  by  the  properties  of  the  integrands  whon  a approaches 
zero  and  when  x approaches  1.  In  thb  first  case,  the  integrand  becomes  infinite  since  !{  be- 
comes infinite  as  nat  log  (//a)  which  necessitate?,  for  numerical  evaluation,  the  transforma- 
tion 3 » tn,  n > 1. 

In  the  following,  the  absolute  value  of  the  function  which  is  defined  by  Equation  [12c] 
or  [13c]  is  denoted  by  F2. 

The  radial  fccce  component  on  the  shroud  is  determined  from  the  axial  velocity  compo- 
nents which  are  induced  at  the  shroud.  Since  this  force  is  normal  to  the  velocity  of  advance 
t?0,  it  does  not  affect  the  transformation  of  power  input  within  the  shrouded  propoilcr  and  is, 
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Figure  10  • Axial  Reaction  Force  Between  Shroud  and  Rotor 


therefore,  not  evaluated.  This  forco  is  significant,  howevor,  for  considerations  of  the 
strength  of  the  shroud. 

3.2  THE  FORCES  ACTING  ON  THE  ROTOR 

As  mentioned  previously,  the  rotor  is  replaced  by  lifting  lines  in  order  to  determine 
the  forces  at  the  rotor.  Then  the  total  force  (including  the  reaction  with  the  shroud)  follows 
from  Kutta's  law,  taking  into  account  the  rosultant  velocity  field  at  a rotor  lifting  line. 

The  tangential  velocity  component  is  generated  by  the  hub  vortex  of  the  rotor  nR  rR, 
and  there  are  no  parts  within  the  tangential  component  either  from  the  shroud  (becauso  ot 
axial  symmetry  of  the  induced  flow)  or  from  the  guide  vanes  (from  the  law  of  Stokes  if  the 
vanes  are  arranged  behind  the  rotor).  Then,  from  Equation  [7],  for  tho  tangential  velocity 
induced  at  the  rotor 


= 


"ll  rH 
4 rir 


>2j3l  ik 

x 2tt 


U51 


With  respect  to  tho  axial  velocity  component  induced  at  a lifting  lino  of  the  rotor,  this 
component  results  from  the  shroud  vortices  and  the  sinks  at  the  disk;  the  latter  includes 
effects  from  tho  lifting  linos  of  the  guide  vanes.  Accordingly 
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(wa)H  = + (u>a)0tS  = ~ +(u»a)/;s 

Within  this  expression,  (u>a)0  D follows  from  Equation  [5]  and  (u>a)D  s from  Figure  2. 
It  is  seon  from  this  diagram  that  (tca)D  s is  a negative  quantity  when  the  sense  of  circula- 
tion at  the  shroud  is  such  that  a pressure  increase  arises  within  the  shroud,  i.e.,  when  y 
is  positive.  Then  (ua)D  0 is  of  the  same  direction  as  vQ  but  (teo)0  s is  opposite. 

With  these  expressions  for  the  induced  velocity  components,  the  force  elements  at 
the  rotor  become 

in  axial  direction  (dT)r  = p(«„/J)[wr  - ( dr  [16] 

in  tangential  direction  ( dQ )H  = p(^/^)[v0  + lf  [171 

Absolute  values  for  the  velocity  components  aro  introduced  into  the  latter  relation  in  order 
to  avoid  mistakes  in  the  signs. 

3.3  THE  FORCES  ACTING  ON  THE  GUIDE  VANES 

Again,  from  Kutta's  law,  these  forces  follow  from  the  velocity  components.  The 
axial  velocity  component  being  continuous,  it  equals  the  axial  component  at  tho  rotor 


The  tangential  force  component  at  the  vanes  which  follows  from  the  axial  velocity 
component  is  without  intorest  for  the  propulsion  of  the  system. 

The  tangential  velocity  component  at  the  vanes  (wt)y  is  determined  from  the  condition 
that  the  guide  vanes  cancel  the  tangential  velocity  field  of  the  rotor  (which,  however,  is 
possible  only  with  an  infinite  number  of  blades.)  Behind  the  rotor,  the  tangential  velocity 
equals  2(w^R,  and  behind  tho  vanes,  it  equals  2{w() y.  Thoreforo,  to  have  the  resultant 
tangential  velocity  zero  in  the  slipstream,  the  following  must  hold 

2 (wt)v  =»  - 2(wt)H 

That  is,  the  total  circulation  at  the  guide  vanes  nyVy  is  equal  but  opposite  to  that  of  the 
rotor  nR  rR,  and  the  circulation  of  the  respective  hub  vortieja  is  equal  but  opposite.  Then, 
the  tangential  component  induced  at  the  vanes  from  the  hub  vortex  of  the  rotor  equals 
2(t c^R  and  that  from  tho  hub  vortex  of  tho  vanes  ; hence,  the  resultant  tangential 

velocity  at  the  vanes  equals 
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l(Wf)yl  = l(wtW  = 


H/tHi  _ nHGn  i’p 
4 rn  x 2 n 


[18] 


From  the  last  relation,  one  obtains  for  an  element  of  the  thrust  at  the  vanes 

( AT)V  - p(nvry)(wt)vdr  [19] 


which  is  in  the  same  direction  ss  the  thrust  element  on  the  rotor. 


U THE  NET  FORCE  OF  THE  SYSTEM  AND  THE  SINK  DENSITY  AT  THE  DISK 

With  the  relations  which  have  been  established  so  far,  expressions  for  both  the  net 
thrust  of  the  system  and  for  the  power  input  can  bo  established. 

The  .set  thrust  T equals 

T = Tft  + 7V  + rs., 

Ts,  from  Equation  [12a],  being  a negative  quantity  for  a positive  y. 

Introducing  Equations  [13c],  [15],  [16],  [18],  and  [19]  pivos  for  the  coefficient  of  the 
net  thrust 


2 n«  _ q */>  p 

w X s v0  2 


[20] 


the  negative  sign  being  necessary  when  introducing  the  absolute  values  \eD\  and  Fa. 

An  expression  for  the  power  input  P is  deduced  from  Equation  [17]  in  connection 
with  Equation  [13c] 


Cg> 


fi/f’-VW 


2 n,tGR  | 
it  A l1  2 


-f-C.s]  [211 


In  order  to  estnblish  an  expression  for  the  sink  density  at  the  disk,  we  refer  to  the 
law  of  momentum  which  states  that  the  time  rate  of  change  of  momentum  equals  the  net 
thrust  of  the  system.  Since  the  scheme  which  has  been  constructed  for  representing  tho 
rotor  flow  leads  to  the  result  that  the  absolute  axial  velocity  equals  |e^|  in  the  final  wake, 
the  change  of  momentum  of  tho  mass  per  unit  time  dm  equals  ( ~Qdm ).  With 

dm  = p2 rndr^v0  + («’<,)/,]  «=  p2rndr[v0  + 1 ^-|  - K«'«)iv<lj 
the  law  of  momentum  yields  the  relation 

T = 2pn £ !<?,,,'[  v0  4 j | - |( u>a)/,,.sl]  rdr 
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In  general,  e r will  be  a function  of  the  radial  coordinate  x.  To  maintain  the  simple 
scheme  of  a circular  sheet  of  vortices  which  is  concentrated  within  the  boundary  of  the  rear 
jet,  it  is  necessary  to  introduce  an  average  of  eD  over  the  radius.  The  reason  for  this  is 
that  a concentrated  circular  vortex  sheet  is  equivalent  to  a uniform  sink  density  over  the 
disk,  as  mentioned  previously.  Otherwise,  to  account  for  a radially  varying  sink  density, 
additional  free  vortex  sheets,  i.e.,  a radially  varying  circulation  at  the  bound  vortices,  must 
be  introduced  which  complicates  the  calculations  considerably. 

The  average  70  is  determined  such  that  the  net  thrust  from  Equation  [20],  which  is 
the  arithmetic  sum  of  the  axial  forces  on  the  substituted  bound  vortices,  equals  the  net 
thrust  from  the  law  of  momentum.  This  leads  to  the  relation 


i£/i  = i/7+  - noli 

1 vu  ' it  A 


[22] 


3.5  RELATIONS  FOR  THE  CIRCULATION  AT  ROTOR, 

SHROUD,  AND  GUIDE  VANES 

When  the  required  net  thrust  of  the  system  is  a given  quantity,  Equation  [20]  repre- 
sents a relation  between  cT  and  the  unknowns  GR  and  Gy  The  same  holds  for  Equation 
[21]  when  the  power  input  is  a given  quantity. 

A condition  for  Gs  is  found  from  the  required  pressure  increase  in  the  plane  of  the 
rotor  Ap  = pR  - pQ.  This  pressure  increase  is  related  to  the  flow  within  which  the  rotor  is 
located  and  whose  axial  component  in  the  plane  of  the  rotor  amounts  to  [t>0  - |(tea)c  s|]. 
They,  it  follows  from  Bernoulli’s  equation  that 


from  which  the  circulation  at  the  shroud  is  represented  by 


123] 


[23b] 


Since  the  induced  velocity  [(wjvj  (l/<?s)]p#s  depends  on  the  radius,  the  required 
pressure  increase  should  be  related  u>  a certain  radius. 

The  loading  coefficient  c7  or  cp  and,  further,  A p/(p/2)v% , nR , and  A are  considered 
given  quantities.  Then  one  obtains  the  circulation  at  the  shroud  Gs  from  Equation  [23b] 
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and,  with  thin  value  for  G s,  the  circulation  GR  from  Equation  [20]  or  [21]  together  with 
[22].  The  equation  for  (nR  GR)  so  obtained  reads  in  terms  of  the  thrust  coefficient 


- («"£*)  * [<Cr-  GsF2)  + + cr  ^ ( cT-?GsFt ) » 0 [24] 

or  in  terms  of  the  power  coefficient 


rGsFi  ( 1 - ^)  + (^)  n*Cp(  1 - GSFJ  - f c>  - 0 125] 


“kG/A 


Gs  F2  \ i nlt  Gk 


The  circulation  at  the  rotor  being  known,  the  circulation  at  the  guide  vanes  equals  the 
total  circulation  at  the  rotor  from  reasons  mentioned  previously.  Therefore 

nvGv  *=  uKGH  [26] 


4.  DESIGN  DATA  FOR  ROTOR,  GUIDE  VANES,  AND  SHROUD 


4.1  ROTOR 

Expressing  the  lift  oo  an  element  of  a rotor  blade  by  the  law  of  Kutta>Joukowsky  and, 
on  the  other  hand,  by  the  lift  coefficient,  and  equating  the  expressions,  the  following  obtains 

(C/.Oft  = 2 & 
rR 

where  the  resultant  relative  velocity  at  the  station  r of  the  blade  is  expressed  by 

VR  “ [ ‘'O  + |^|  ~ l(W«>Orsl]  + K “ 

Nondimensionally,  these  expressions  read 


The  density  eD/uQ  is  known  from  Equation  [22],  the  axial  induction  from  the  shroud 
|[(to0/o0)  (l/ff5)]0tS  | from  Figure  2. 

In  connection  with  the  last  two  relations,  the  direction  of  the  relative  velocity  by 
which  the  position  of  the  section  is  fixed  follows  from  tho  relation  (see  Figure  11) : 
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Figure  11  • Velocity  and  Force  Diagram  at  a Section  of  the  Rotor 


Further,  in  the  design  of  the  'otor,  its  local  cavitation  number  must  be  known.  This 
is  given  by  the  expression 


CONFIDENTIAL 


27 


CONFIDENTIAL 


where  aQ  is  the  cavitation  numbor  of  tho  undisturbed  flow. 

Cavitation  will  not  toko  place  on  a section  of  tho  rotor  if 


Ok  § 


P~Po 


£.  tv 
2 0 


min 


whore  tho  term  on  tho  righUhand  side  denotes  tho  tbsolute  valuo  of  tho  nondimonsional  mini* 
mum  of  tho  pressure  on  tho  section. 

By  moans  of  these  formulas,  the  design  of  the  rotor  follows  usual  procedure,  viz.,  to 
determine  c L and  l so  that  the  calculated  product  (c^  t)  is  satisfied  and  that,  i i the  same 
time,  the  onset  of  cavitation  is  avoided.  This  roquiros  knowledge  of  both  lift  versus  angle 
of  attack  curves  and  pressure  distribution  curves  of  sections  in  cascado.  On  the  other  hand, 
the  knowledge  of  two-dimonsional  cascado  effocts  on  the  sections  is  sufficient  for  this  prob- 
lem  since  the  flow  at  the  rotor  for  a constant  circulation  is  essentially  two-dimensional  (in 
contrast  to  the  unshroudod  propolicr  whoro  the  flow  is  ossentiallx  three-dimensional). 

4.2  GUIDE  VANES 

Analogously  as  for  tho  rotor,  there  is  obtainod 


i.  I \ _ <»  fy  _ "h  T); 

(C,/)V  - l yv  ~ i-  y- 


since  the  condition  of  Equation  [26]  must  be  true  in  order  to  cancel  tho  tangential  velocity 
field  of  tho  rotor  by  means  of  tho  vanes,  in  this  flow  condition,  tho  resultant  relative  veloc- 
ity at  tho  vanes  Vf,  is  expressed  by 

-[*+!*!-  + <«v>* 

Nondimensionally,  the  last  two  relations  become 

"N  I Vo  \ /; 


The  angle  ftv  of  tho  resultant  relative  velocity  at  the  vanes  with  the  plane  of  tho 
rotor  is  determined  from  the  following  equation  (see  Figure  12) 


tun  = 2 rrx 


. + !&'_/•  |/r* -U  I 

■i  i-.  '*11  v„ 


nR  °R 


[31] 


rtfttirinrtjTi  is 
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Figure  12  • Velocity  and  Force  Diagram  at  a Section  of  the  Guide  Vanes 
For  die  local  cavitation  number  at  the  vanes,  the  following  obtains 


if  the  approximation  pv  - pR  is  applied. 

4.3  SHROUD 

In  accordance  with  thin  airfoil  theory,  the  shroud  section  is  considered  to  consist 
of  an  infinitely  thin, camber  line  section  (which  produces  circulation)  and  of  a superimposed 
thickness  form,  the  combined  effect  of  which  on  the  velocity  distribution  is  obtained  when 
adding  the  individual  effects.  In  the  following  treatment,  the  shape  of  the  camber  line  and 
its  position  relative  to  the  undisturbed  flow,  i.e.,  its  geometric  angle  of  attack,  are  first  de- 
termined so  that  an  elliptic  distribution  of  circulation  is  realized  on  the  shroud,  and  then  a 
thickness  function  in  axisymmetrical  flow  is  considered. 

The  shape  of  the  camber  line  follows  from  the  boundary  condition  that  the  velocity 
vector  at  the  camber  line  is  tangent  to  this  line,  i.e.,  that  this  line  represents  a streamline 
of  the  flow.  The  flow  at  the  camber  lino  is  a result  of  the  velocity  of  approach,  the  self- 
induction  of  the  shroud,  and  the  inductions  from  both  the  sink  disk  and  the  sheft  sink.  Con- 
sequently, the  axial  component  oi  the  resultant  flow  amounts  to 

*o  + («'a)i.rS  + ( u>a)SpD  + <tt-o).SffA 
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and  the  radial  component  to 


< "v\S(.s  + < «’r).s,/. 


The  camber  lino  , with  coordinates  a and  (r  - U),  is  a streamline  if 


dr  _ ( »<y)s.,s  + (ti  vjs»/i  tjA 

r//7  ~ C0  + + ( «*,)>•,* 


In  this  relation,  the  induction  from  the  shaft  sink  has  boon  allowed  for  in  order  to  have  tho 
continuity  oquntion  satisfied. 

introducing  tho  nondimensional  coordinates 


U 


r_-  A* 

Hi ' 


and 


V/a 


gives  a first-order  differential  equation  for  tho  shape  of  tho  camber  line  y » y(z) 


dy  = 


I + 


kL*  4 'fK  gl 


dr 


[33] 


in  which  relation  the  terms  aro  determined  by  tho  equations  or  diagrams  as  listod  in  the  table 
on  pago  30.  The  tablo  also  indicates  the  sign  of  tho  respective  volocity  components  for  a 
positive  circulation  at  the  shroud  and  for  a sink  density  at  tho  disk. 

From  Equations  [3]  and  [4],  the  integrands  of  both  («aa)s,p  and  (u,f)s  p become  in- 
finite for  x = 1 when  z becomes  zero,  i.e.,  when  the  point  of  reference  coincides  with  the 
edge  of  the  sink  disk.  These  infinities  ariso  from  the  assumption  that  the  vortices  of  the 
shroud  are  situated  on  the  disk  cylinder  instead  of  being  arranged  on  tho  camber  line.  This 
means  that  in  the  neighborhood  of  z » 0,  the  finito  distance  d/2  betwoen  cambor  lino  and 
disk  cylinder  should  be  taken  into  account  although  it  can  be  noglected  outside  of  a certain 
interval  around  z = 0.  For  this  roason,  the  induction  at  tho  shroud  from  tho  disk  as  repre- 
sented on  Figures  7 and  8 can  be  considered  accurate  in  an  interval  of  z botween  - 1 and 
z about  - 0.2,  the  last  figure  being  obtained  from  comparative  calculations  of  the  velocity 
components  at  lines  r ■ K + d/2  and  r = K. 

Bbvwsen  z - +0.2  and  z = -0.2,  tho  induced  velocity  components  as  represented  on 
Figures  7 and  8 are  exaggerated.  With  respect  to  the  integrand  of  Equation  [33]  which 
must  be  known  within  tho  entire  interval  from  s = -ltoa«=+l,u  cloco  approximation  is 
obtained  when  reading  off  (irf)s  /;  botweon  z = -0.2  and  z - 0 on  tho  dotted  curves  of  Fig- 
ure 8.  For  these  dotted  curves,  tho  radial  distance  of  the  points  of  reforen.' . equals 
(iif  + d/2),  i.e.,  * = 1 + ( d/l)h ; those  poii  t r.  are  vory  close  to  tho  cambor  lino  around  z = 0. 
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To  silow  for  the  influence  of  the  thickness  of  the  shroud,  the  calculations  have  been  made 
for  d/l  - 0.04  and  d/l  ■ 0.06  (Figure  8a);  these  curves  are  so  related  that  equal  coordinates 
are  obtained  for  equnl  quantities  (d/l)  h. 

For  (wJs'j),  this  component  equals  zero  at  the  point  with  coordinates  a ~ 0 and 
x • 1 + (d/l)h.  Since  this  component  occurs  within  a sum,  together  with  the  relatively  great* 
er  number  1,  it  is  sufficiently  accurate  to  interpolate  this  sum  between  a - -0.2  and  a - +0.2 
from  the  known  values  at  a - -0.2  and  o ■ +0.2. 

In  this  way,  a closer  approximation  to  the  velocity  Held  at  the  camber  line  of  the 
shroud  is  obtained  than  by  replacing  the  rink  disk  by  a point  disk  at  the  axis,  the  strength 
of  which  is  such  that  the  sink  disk  and  the  point  sink  are  acted  on  by  an  equal  interaction 
force  from  the  shroud.  This  approximation  has  been  proposed  by  Horn  for  the  purpose  of  de- 
termining the  circulation  at  a Kort  nozzle. s Tho  two  velocity  fields,  however,  differ  groatly 
in  the  vicinity  of  the  shroud  and  becomo  sufficiently  coincident  only  at  a great  distance  from 
the  disk. 

The  integration  of  Equation  [33]  requires  the  determination  of  an  integration  constant 
from  tho  given  ordinate  of  the  camber  line  at  an  arbitrary  station  a.  Corresponding  to  the 
geometric  configuration  of  tho  propulsion  system,  this  constant  is  chosen  so  that  y - d/2 
for  a « 0.  With  this,  a first  approximation  for  tho  shape  and  for  the  geometric  angle  of  attack 
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of  the  camber  lino  is  obtoinod  when  integrating  Equation  [33]  by  numerical  mothods.  A soc- 
ond  approximation  fellows  when  arranging  tlio  vortices  of  the  shroud  on  tho  first  approxima- 
tion for  the  combor  lino  and  whon  calculating  tho  velocity  components  which  occur  in  Equa- 
tion [331  on  the  cambor  line  instoad  of  on  tho  rotor  cylindor.  The  previously  mentioned 
tables  by  Kueehemann  aro  very  useful  for  this  purposo.*  Tho  velocity  components  which 
are  induced  at  the  first  approximation  of  the  shroud,  both  from  tho  sink  disk  and  from  the 
shroud  itself,  can  bo  ascertained  within  a reasonable  amount  of  time  by  applying  these  tables. 

The  maximum  velocity  (or  the  minimum  pressure)  on  the  shroud  must  bo  known  in  order 
to  determine  the  cavitation  number  of  tho  shroud.  Following  tho  procedure  with  thin  airfoils, 
a close  approximation  to  the  velocity  distribution  on  the  shroud  is  obtainod  when  adding  the 
velocities  on  the  infinitely  thin  camber  line  section  and  those  on  the  thickness  shape.  An 
additional  velocity  distribution  associated  with  angle  of  attack  need  not  bo  considered  sinco 
the  shroud  is  supposed  to  work  in  the  shock-free  flow  condition. 

By  definition,  the  thickness  form  produces  no  lift  at  r.oro  angle  of  attack.  Thorefore, 
in  two-dimensional  flow,  it  is  symmetrical  about  the  direction  of  tho  undisturbed  velocity. 

For  tho  same  reason,  tho  velocity  field  of  the  thicknoss  shopo  can  be  represented  by  a suit- 
able sink-source  distribution  on  tho  axis  of  symmetry  of  tho  section.  This  mothod  is  also 
applicable  to  the  thicknoss  function  of  a ring-shaped  hydrofoil  (in  contrast  to  tho  mothod  of 
conformal  mapping  which  is  restricted  to  two-dimensional  flow).  Kueehemann  has  investi- 
gated the  effoctof  axial  symmetry  in  this  way,  assuming  ring-shopod  sink-source  distribu- 
tions by  which,  in  the  case  of  two-dimensional  flow,  symmetrica!  Joukowsky  profiles  ore 
generated.6  In  the  axi symmetrical  case,  the  generated  annular  sections  are  no  longer  sym- 
metrical about  the  direction  of  tho  undisturbed  flow  but  havo  a curved  middle  lino;  the  cambor 
increases  when  the  ratio  of  thickness  to  chord- length  d/l  increases.  Tho  influence  on  the 
pressure  distribution  is  shown  by  Kueehemann  on  annular  half-bodies  (which  are  generated 
by  a ring  source  within  the  undisturbed  flow)  from  which  it  follows  that  the  prossure  distri- 
bution on  the  body  is  no  longer  symmetrical  in  a meridian  piano  but  tho  peak  of  suction  is 
greater  inside  than  outside,  tho  asymmetry  increasing  with  thicknoss. 

It  follows  from  those  investigations  that  whon  taking  a thickness  form  whose  pressure 
distribution  is  known  in  two-dimensional  flow  (e.g.,  one  of  the  NACA  basic  thickness  forms), 
we  cannot  rely  immediately  on  tho  two-dimensional  pressure  distribution  ia  sxi symmetrical 
flow.  In  order  to  ascertain  the  difference  of  tho  pressure  distribution  for  the  same  section  in 
these  two  cases,  it  has  been  attempted,  by  moans  of  Equations  [3]  and  [4],  to  determine  a 
distribution  of  ring-shaped  sinks  end  sources  on  the  axis  so  that  the  thickness  form  becomes 
a streamline.  In  the  analogous  twxt-dimensionai  problem,  the  sink-source  distribution  follows 
from  an  integral  equation.  In  the  case  of  axial  symmotry,  it  has  not  been  possible  yet  to 
establish  the  corresponding  integral  equation  bccauso  of  complications  wnich  arise  from  tho 
expressions  [3]  and  [4], 

Although  present  knowledge  does  not  pormit  comparison  of  sections  of  equai  shape 
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in  the  two  flow  conditions,  an  idea  of  the  order  of  magnit.do  of  tho  influence  of  axial  sym- 
metry is  obtained  when  comparing  sections  of  equal  sink*sourco  distribution.  Since  the  in- 
fluence on  both  the  shape  and  the  prossuro  distribution  increases  with  increasing  thickness, 
it  is  of  interest  to  know  whether  or  not  a thickness  can  bo  determined  below  which  these  in- 
fluences are  so  smaii  that  they  can  bo  neglected.  Taking  tho  afore-mentioned  sink-source 
distributions  wlJch,  in  the  two-dimensional  case,  lead  to  symmetrical  Joukowsky  profiles, 
the  results  are  as  follows:  Bolow  d/l  about  0.08,  the  differences  in  pressure  distribution  as 
compared  with  the  two-dimensional  section  of  equal  thickness  chord-length  ratio  are  negli- 
gibly small  in  the  range  of  h investigated  (0.7  to  2.1).  Further,  within  this  range  of  d/i  and 
A,  the  deformation  of  the  contour  is  allowed  fo*  with  sufficient  act  jracy  by  a curvaturo  of  the 
middle  tine  of  the  section  without  appreciable  change  of  the  ordinates,  the  latter  being  re- 
ferred to  the  middle  line.  The  curvature  is  small  and  in  the  interval  investigated,  it  depends 
linearly  on  d/l,  increasing  slightly  when  A increases  (Figuro  13). 


<t 

T 


Figure  13  • Camber  of  Thickness-Forms  in  Axi symmetrical  Flow 

It  is  concluded  from  these  calculations  that  the  two-dimensional  pressure  distribution 
of  a thin  thickness  form  (d/l  being  not  greater  than  about  0.08)  can  bo  applied  for  axisymme- 
trical  flow  if  the  thickness  form  is  given  a slight  curvature  corresponding  to  Figure  13. 

An  expression  for  the  minimum  pressure  on  the  shroud  can  be  deduced  on  this  basis. 
Within  the  assumptions  of  linearized  thin  airfoil  theory,  the  velocity  at  any  point  of  the 
shroud  is  the  sum  of  the  veiocities  from  the  circulation  distribution  on  the  camber  line  and 
from  the  thickness  form:  V * V s*  V T. 

The  latter  is  assumed  to  bo  known,  e.g.,  from  one  of  the  NACA  thickness  forms;  the 
former  is  the  velocity  just  outside  or  inside  of  the  camber  line,  i.e.,  just  outside  or  inside 
of  tho  vortex  sheet.  Let  u>0  be  tho  velocity  at  a point  of  the  outside  surface  with  axial 
coordinate  a,  and  t ci  the  velocity  at  a point  of  the  inside  surface  with  the  same  axial  coor- 
dinate. Then 
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tv0  ~ W;  = y 


where  y is  the  circulation  per  unit  length  at  a.  Further,  just  on  the  outside  of  the  vortex 
cylinder,  the  relation 


M*  . 

Sheet 


+ «•' 


holds,  and  just  on  the  inside  at  a point  of  equal  a 


uv 


u> 


Sheet 


- w' 


where  to' is  the  increment  of  velocity  as  compared  to  the  velocity  within  the  vortex  sheet  at 
equal  a.  From  those  relations,  it  follows  for  the  unknown  velocity  to' 


w‘  = 


Y_ 

2 


Then,  at  a pc'  t on  the  outside  of  the  camber  line  where,  with  & positive  y,  the  smallest 
pressures  •bund,  the  following  obtains 


The  angle  a is  the  declination  of  the  camber  tine  with  respect  to  the  z*oxis  which  is  small. 

Relative  to  the  values  of  D and  of  («r)s  D in  the  interval  of  a between  a = -0.2 
and  s • +0.2,  see  the  explanation  given  in  connection  with  Equation  [33], 

The  velocity  distribution  VT/v0  is  known  for  the  thickness  form.  Then  the  maximum 
of  the  velocity  at  a shroud  of  finite  but  small  thickness  becomes 

M -(*  + &) 

\ V0  Va‘mut 


In  order  to  avoid  the  onset  of  cavitation  on  the  shroud,  the  following  relation  must  be  satis* 
fled  for  the  cavitation  number  of  the  oncoming  flow 


g 


1 
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5.  INFLUENCE  OF  VISCOUS 

From  Figure  11,  it  follows 
• denoted  by  a prime 

( dT V - ( dT)R 
= (dT)H  - «, 

UQ),;  “ (dQ)H  + (dD)^  cos  0H  = (dQ)K  (l  + ~^) 

Within  these  expressions,  is  the  drag-lift  ratio  at  the  section  of  the  rotor  under  considera- 
tion. 

Analogously,  one  obtains  from  Figure  12  for  a section  of  the  guide  vanes 

MTV  = UT)V  - (dD)y  sin  $v  = ( dT)v  ( 1 - e tan  0V ) 

At  tho  shroud,  the  viscous  drag  of  the  shroud  and  the  (nonviscous)  interaction  forces  are  in 
the  same  direction  when  y is  a positive  quantity.  In  this  case,  the  total  drag  of  the  shroud 
becomes 

Ts  *=*  Ts  + Ds 

or,  nondimensionally, 


FLOW  ON  THE  FORCES  AND  ON  THE  EFFICIENCY 

for  the  rotor  when  the  respective  forces  in  viscous  flow  are 

- ( dD)H  sin  = (dT)u  - tn(dL)H  ain 


(dT), 


ain0*  - (dT)H  ( 1 - e*tan  0*) 


{ Ct)s 


Ts 

%(R*-r£)Trv0t 


(Cjr)s  + ( CD) g 


4 h 

1 -V 


if  (cD) s,  as  usual,  is  referred  to  the  surface  of  the  shroud  2 Rnl. 

The  forces  on  rotor  and  vanes  can  be  integrated  when  it  is  assumed  that  the  drag-lift 
coefficients  are  independent  of  tho  radius.  Introducing  relations  [16],  [17].  and  [19]  for  the 
nonviscous  parts  of  the  forces,  and  relations  [28]  and  [31]  for  the  angles,  one  obtains 


<«,>;  - w,  - <4  [( ■ * l)< 1 - - c»0  t,  iU  H w 


(c,)'  - (V  + *4  i iSL  [i , 1 - V)  - ^ ( 1 - *,)] 


Gk 
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«r>v  * <«rV  - + ' - V-  O,0S  A)  J<0  [371 

It  follows  for  tho  coefficient  of  the  net  thrust  c j that 

ct  = ( crV  + ( crV  “ < crV  = ct~(  £k  + (\)4„  i '*3ft 


[( 1 + -2  !$!)'  >-»>-  *0®  £),„>]  - «•*  r^h 


138] 


With  these  expressions  for  cj and  Cp , tho  efficiency  of  the  system  becomes 

CT 

*1  ^ «f. 

where  the  ideal  efficiency  rji  is  obtained  when  the  expressions  [20]  and  [21]  for  cT  and  cp 
are  introduced 
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[39] 


For  Gs  = 0 (unshrouded  propeller),  this  expression  for  jj.  equals  that  from  simple  momentum 
theory.  This  is  necessary  because  of  the  assumptions  that  the  guide  vanes  cancel  the  tan- 
gential volocity  field  entirely  and  that  the  circulation  is  independent  of  the  radius;  these 
are  just  the  assumptions  of  simple  momeruim  theory. 

For  the  reduction  of  the  ideal  efficiency  due  to  drag,  one  obtains 


Assuming  tho  conditions  of  an  unshrouded  propeller  without  guide  vanes,  vis., 

Gs  = 0,  (0^)5  ■ 0,  « y = 0 and,  further,  putting  xh  = 0,  it  follows  from  the  last  relation  that 
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The  third  term  in  the  denominator  is  small  compared  to  1 and  can  bo  neglected. 
Further, 


Then,  from  the  above  expression,  the  well-known  approximate  formula  r t)(  of  an  unshrouded 
propeller  which  has  boen  deduced  by  von  Kdrmdn  and  Bienen  is  obtain*  *1 : 


n = 1 fLfitAi 

" -It 


It  is  difficult  to  draw  genet*.  .sions  from  the  relations  for  the  efficiencies.  For 
a special  case,  viz.,  that  for  the  low-speed  pumpjet  for  DD710,  Mrdel  3246,  the  following 
design  data  are  given 


e'H  = 1.08,  A = 0.45,  = 0.3,  h = 1.2 

r*  2 
2 l’° 

By  means  of  successive  approximations  as  described  in  section  7 of  this  report,  the  follow- 
ing obtains 


- (fl ” °-022 

For  determining  those  drag-lift  coefficients,  the  drag  coefficients  (cD)R  - 0.010  and 
(cD)v  = 0.024  have  been  introduced,  corresponding  to  the  respective  Reynolds  numbers  and 
angles  of  attack.  The  latter  figure  is  greater  th^sn  the  former  since  the  sections  of  the  vanes 
work  at  a great  lift  coefficient  with  which  a groat  pressure  resistance  is  connected.  In  addi- 
tion, the  Reynolds  number  of  the  sections  of  the  yanes  is  smaller  than  that  of  the  rotor;  be- 
cause of  this,  the  frictional  coefficient  becomes  greater. 

Assuming  (cfl)s  - 0.013  (which  figure,  however,  is  uncertain)  and  using  Equations 
[39]  and  [40],  the  following  values  are  derived 


t)  j = 0.808,  = 0,841,  q » 0.679 

From  model  tests,  i;  has  been  determined  as  0.656.  Further,  a value  of  0.77  is  obtained 
from  these  calculations  for  the  ratio  of  the  net  thrust  to  the  thrust  at  the  rotor  and  a value 
of  0.725  from  the  model  tests.  These  differences  between  calculations  and  test  results  may 
arise  from  the  uncertainty  of  (cfl )s  which  probably  is  greater  than  assumed  becauso  of  sep- 
aration which  arises  from  the  shape  of  the  shroud  and  from  the  pressure  increase  at  the 
rote r. 
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6.  INFLUENCE  OF  A FINITE  NUMBEL  OF  BLADES 


The  schomo  which  has  been  introduced  for  the  action  of  tho  rotor  has  Leon  chosen  as 
simple  as  possible,  viz.,  constant  circulation  along  the  lifting  linos,  a hub  vortex  of  the 
combined  strength  of  the  lifting  lines,  and  a system  of  vortices  within  the  boundary  of  tho 
rear  jet.  Furthermore,  it  has  been  assumed  that  the  number  of  blades  is  infinite.  In  this 
section,  the  last  assumption  will  not  bo  made  but  the  velocity  components  which  are  induced 
at  a lifting  line  by  nR  helical  tip  vortices  togother  with  a straight  hub  vortex  of  strength 
(ftJl  f ft)  be  considered. 

Without  going  into  details  of  the  deduction  (e.gM  seo  Kawada7),  it  is  stated  that  the 
ratio  of  both  the  tangential  and  axial  velocity  components  at  a lifting  line  for  a finite  number 
of  blades  to  the  respective  quantity  for  on  infinite  number  of  blades  is  expressed  by  the 
following  relation 


whero 


w,  nv,  j 2 nK  s 

(wt)  (w,j  tan  0R 


s ~ 2.  K* 


( mnk  \ 
\ tan  I 


Within  these  expressions,  tan  (JR  is  the  pittu  ou^lu  of  the  tip  vortices,  / and  K are  the 
modified  Bessel  functions  of  the  first  and  sect  r.d  kind,  respectively,  and  the  prime  indicates 
the  derivative  with  rospect  to  the  argument  W.e..  \arying  tho  number  of  blades,  it  is  assumed 
that  both  the  total  circulation  (nR  VR)  and  tho  pitch  anglo  flR  are  kept  constant.  This  im- 
plies that  tho  loading  coefficient  varies  when  the  numbor  oi  blades  is  varied. 

The  sum  S depends  on  the  pitch  angle  pR  and  on  tho  number  of  blades  nR.  For 
nR  - 5,  the  quantity  i(  which  is  (olatca  to  S by 


i,  1 trfU  fj.. 
2 x - i 


is  represented  on  Figure  14.  This  quantity  t,  is  denoted  as  tho  “induction  factor.”  Then, 
the  ratio  of  the  velocity  i.ompononts  becomes  in  tonus  of  tho  induction  factor 

_L?  . , | lr  £ 

\&t>  iu-n)  1 ‘ 1 

Bj  the  use  of  singlo  up  »'ortice  , tie  mil.  *t«l  • ».lm  tty  components  are  exaggerated 
near  the  tip  since,  when  * - I,  both  u.,  uud  Loconic  infinite.  Beyond  a certain  distance 
from  the  tip,  however,  this  formula  is  suitable  tor  obtaining  an  order  of  magnitude  for  the 
influence  of  a finite  numbci  blades.  For  / ih  .53  liogtoc.-,  c.g.,  tho  right-hand  side  bo- 

cornos  l.u^O  at  x - o.C,  ■>!..  ...  ..  ii  ..  S.-.1  in f utung  in  practice,  tho 

influence  .f  a finite  iiumhet  ( f 1 lade.- i • i.f  >.  .'ey  i.i.j.i.i Thir>  holdt  for  the 
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fl,t  m digrtis 


Figure  14 


assumption  that  the  circulation  at  the  lifting  lines  is  independent  of  the  radius  but  is  differ* 
ent  when  the  circulation  distribution  depends  on  the  number  of  blades  (ns  it  does  with  an  un* 
shrouded  optimum  propeller). 

The  same  considerations  are  applicable  for  guide  vanes  with  a finite  number  of  blades. 
The  influence  of  the  number  of  blades  on  the  induced  velocity  component  is  greater  for  the 
guide  vanes  than  for  the  rotor  since  fiv  > fiR . The  induction  factor  increases  when  jS  increas- 
es, hence  the  factor  is  somewhat  greater  for  the  guide  vanes  than  for  the  rotor.  For  /3W  - 33 
degrees  and  \ *>  0.45,  e.g.,  $v  becomes  53  degrees,  with  which  quantity  the  correction  to  the 
velocity  components  from  an  infinite  number  of  blades  is  1.028  for  tiy  •*  6 and  x *»  0.6.  When 
the  number  of  blades  is  increased,  the  percent  correction  decreases  fairly  rapidly,  roughly 
to  half  its  value  when  one  blade  is  added  (o.g.,  to  1.016  for  nv  « 7). 
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7.  METHOD  OF  DESIGN 


The  method  of  design  of  a shrouded  unit  from  the  foregoing  considerations  is  analo- 
gous  to  that  of  the  application  of  the  usu?J  propeller  theory  for  design  purposes. 

The  given  quantities  are  loading  coefficient  e£  (or  c£),  advance  coefficient  A,  pressure 
increase  at  the  rotor  p R/pQ , and  the  cavitation  number  of  the  oncoming  flow 
ao  = ( Po  ~ P*)  / £ vf  When  a value  of  A - 1/2R  is  chosen,  the  circulation  at  the  shroud, 
required,  in  order  to  realize  the  pressure  increase  oecomes  known  from  Equation  [23b]  and 
Figure  2. 


[23b] 


The  quantity  Ap/(p/2)?-j  is  related  to  pR/p0  and  aQ  by 


The  next  step  is  to  determine  the  loading  coefficient  of  the  system  when  operating  in 
nonviscous  flow;  this  requires  the  determination  of  the  drag-lift  coefficients  of  the  compo- 
nents of  the  system.  This  problem  can  be  solved  by  successive  approximations  when  putting 
in  a first  step,  Cj  **  cT  (with  the  power  coefficient  given,  the  method  is  entirely  analogous). 
Then,  the  circulation  at  the  rotor  follows  from  Equation  [24] 


(-^f  " n(er-  GsF ’ + + c*  ?*<cr  “ 2 GsFt)  = 0. 


F2  is  known  from  Figure  10  as  a function  of  A. 

Putting  Cj  - c~,  a first  approximation  is  obtained  for  (nR  GR)  with  which  a first 
approximation  for  the  sink  density  follows 


+ 2 
7T  A 


1 


[22] 


At  this  point, « check  is  mode  as  to  whether  or  not  the  quantity  A chosen  together 
with  Uie  given  quantities  are  consistent  with  a cavitation-free  flow  on  the  shroud.  The 
condition  for  this  is  that 


- 

' VqIimi 
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where 


(Relative  to  (tca)s  q,  see  the  discussion  on  page  17  for  positive  values  of  a.) 

In  Equation  [34a],  the  cavitation  number  of  the  oncoming  flow  oQ  is  a given  quantity. 
Further,  Vf/vQ  is  known  if  a certain  thickness  form  of  the  shroud  is  selected.  Within  Equa* 
tion  [34],  Gs  and  |e0|  have  already  been  determined  and  the  velocity  functions  are  known; 
see  the  table  on  page  30.  In  these  first  steps,  the  declination  of  the  camber  lino  a is 
assume!  zero. 

If  the  value  of  A selected  does  not  satisfy  Equation  [34a],  a greater  length  of  the 
shroud,  i.o.,  a greater  A,  is  chosen  and  the  calculations  are  repeated. 

Assuming  that  the  A satisfies  the  cavitation  condition  of  Equation  [34a],  first  approxi- 
mations of  the  cavitation  number  on  both  the  rotor  and  the  guide  vanes  are  obtained  from 


and  from 


where 


and  where 


[29] 


[32] 
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In  these  first  approximations,  both  a ^ and  uv  aro  calculated  oni\  ;•!  il'i  • ..»! 1 1 - 
x » 0.7,  which  is  sufficient  for  determining  the  drag-lift  coefficients. 

The  maximum  values  of  tho  lift  coefficients  which  aro  permissible  from  the  point  of 
view  of  onset  of  cavitaticn  on  both  tho  rotor  and  tho  guide  vanos  follow  from  these  cavitat;on 
numbers  together  with  diagrams  on  tho  critical  capitation  number  or  on  tho  minimum  pressure 
of  a family  of  suitable  sections.  In  ordep  a obtain  approximations  for  tho  drag-lift  coefficients, 
the  drag  coefficients  (cD)K  and  ( cD)v  must  bo  known  for  tho  typo  of  sections  chosen.  A 
reasonable  first  assumption  for  tho  rotor  is  cp  = 0.008.  For  the  vanos,  cD  is  obtained  as 
the  sum  of  the  frictional  coefficient  and  the  prossuro  drag  coefficient,  the  latter  doponding  on 
the  angle  of  .attack  which  is  usually  great  at  the  vanes.  At  the  shroud,  the  drag  coefficient 
(cD)s  may  be  groater  than  the  frictional  coefficient,  in  spito  of  the  nhock-freo  flow,  as  a 
consequence  of  separation  inside  of  the  shroud  which  arises  from  its  action  as  a diffusot 
and  from  the  action  of  the  rotor.  No  information  on  this  effect  could  be  found  in  the  literature. 
Measurements  on  tho  drag  coefficient  of  annular-shapod  wings  aro  restricted  to  acceloratod 
flow,  i.e.,  the  nozzles,  in  this  case  and  for  a thickness  ratio  of  the  section  of  about  20  per- 
cent, a value  of  0.015  for  cD  at  Reynolds  number  4.10  s has  been  determined.  This  figure  is 
somewhat  greater  than  would  be  expected  for  a two-dimensional  wing  or  equal  section. 

The  drag-lift  coefficients  of  tho  components  of  the  system  being  known  approximately, 
the  loading  coefficient  of  the  system  in  nonviscous  flow  cT  follows  from 


The  approximation  is  repeated  in  a second  step  with  this  value  of  cT.  Subsequent 
steps  aro  necessary  until  the  value  of  cT  which  results  in  a certain  step  does  not  differ 
appreciably  from  that  which  was  assumed  for  that  step. 

When  cT  has  been  determined  in  this  way,  tho  exact  values  of  Gs,  Gr,  tD/v0, 
(Vr/v 0)  and  (Vy/v0) can  be  ascertained,  the  last  two  quantities  as  functions  of  tho  radius 
x.  Then  the  products  (c^  t)  at  both  rotor  and  vanes  follow  as  functions  of  a:  from 


and  from 


[30] 
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Further,  at  each  radius,  tho  cavitation  numbers  aR  and  av  become  known  from  Equa- 
tions [29]  and  [32),  respectively.  With  the  aid  of  tho  afore-montionod  diagrams  for  critical 
cavitation  numbers  of  families  of  sections  (which  should  include  the  two-dimensional  cas- 
cade effect),  the  products  (cl  l)  are  so  split  up  into  their  factors  cL  and  l that  c L equals  or 
becomes  smaller  than  that  lift  coefficient  for  the  respective  local  cavitation  number  which 
is  permissible  relative  to  the  onset  of  cavitation.  With  this  lift  coefficient,  the  angle  of 
attack  of  tho  respective  section  against  tho  resultant  relative  velocity  is  determined  when 
the  lift  versus  angle  of  attack  curves  of  the  sections,  including  the  c&scado  effect,  aro 
known.  The  direction  of  the  resultant  relative  velocity  against  the  plane  of  the  rotor  is  ob- 
tained for  the  rotor  from 


tan  = 


£ _ "fe-ff/t 
k 2rrx 


128] 


and  for  the  guide  vanes  from 


tXTifty  ~ 


2jtX 
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The  design  of  the  shroud  requires  tho  integration  of  Equation  [33]  for  the  determina- 
tion of  the  shape  of  the  camber  line  and  its  orientation  with  respoct  to  the  axis.  In  this 
equation,  all  quantities  ere  known  from  the  preceding  numerical  calculations.  Tho  influence 
of  axi symmetrical  flow  on  the  pressuro  distribution  of  the  thickness  form  is  approximately 
compensated  for  by  an  additional  camber,  from  Figuro  13,  to  be  superimposed  on  the  shape 
which  follows  from  an  integration  of  Equation  [33].  At  least  two  steps  are  necessary  for 
determining  the  shape  and  the  geometric  anglo  of  attack  of  the  camber  lino  from  Equation 
[33],  viz.,  arranging  tho  vortices  first  on  tho  rotor  cylinder  and,  afterwards,  on  the  camber 
line.  Tables  by  Kuechemann  in  which  ti.o  velocity  components  from  both  vortex  and  sink 
rings  have  been  tabulated  are  very  useful  for  the  second  and  higher  approximations. 8 

Finally,  the  efficiency  of  the  system  is  obtained  from  Equations  [39J  and  [401.  In  the 
latter  relation,  the  drag-lift  ratios  <R  and  tv  aro  considered  independent  of  r.  In  general, 
it  will  be  sufficient  to  introduce  the  respective  quantities  at  z « 0.7  as  a suitablo  average. 


8.  CONCLUSION 

The  considerations  of  this  p-ipor  are  based  on  a circulation  at  the  bound  vortices  of 
the  rotor  which  is  independent  of  the  radius.  This  case  represents  the  optimum  with  respoct 
to  efficiency  for  a shrouded  propeller.  Tho  flow  for  circulation  distributions  which  differ 
from  the  optimum  can  bo  determined  in  principle  from  tho  offocts  of  tho  free  vortex  shouts 
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which  are  a consequence  of  a radially  varying  circulation  and  from  their  interference  with 
the  shroud  vortices.  Such  problems  arise  when,  e.g.,  a shrouded  unit  is  given  which  satis- 
fies the  optimum  circulation  in  the  design  condition  and  when  the  circulation  distributions 
of  both  rotor  and  shroud  are  to  bo  determined  in  an  off-design  condition.  For  a finite  number 
of  blades,  this  problem  leads  to  an  integral  equation  which  can  only  be  solved  by  approxi- 
mate methods.  The  velocity  fiela  of  vortex  sheets  with  an  arbitrary  pitch  distribution  must 
be  known  in  order  to  obtain  a solution.  The  necessary  numerical  quantities  of  this  velocity 
field  have  been  ascertained  at  the  lifting  lines  for  a range  of  blade  numbers.  The  determina- 
tion of  the  induced  velocities  at  the  shroud,  however,  would  require  a great  deal  of  addition- 
al numerical  work 

In  addition  to  this  theoretical  work  on  units  with  a radially  varying  circulation,  it  is 
the  opinion  of  the  author  that  the  following  experimental  work  is  necessary  for  a further  de- 
velopment of  shrouded  propellers: 

a.  Studies  of  pressure  distribution,  separation,  and  radial  distribution  of  the  inside  flow 
of  annular- shaped  wings  whose  flow  is  retarued.  The  available  experimental  papers9,10 
give  information  on  pressure  distribution,  but  no  systematic  information  could  be  found  either 
on  the  drag  or  on  the  velocity  distribution  of  the  inside  flow.  In  addition,  the  influence  of 
the  shaft  on  the  afore-mentioned  quantities  should  be  ascertained. 

b.  Studios  on  shrouded  propellers  with  a transparent  shroud  are  needed  in  order  to  in- 
vestigate the  cavitation  performance  of  shroud,  rotor,  and  guide  vanes.  Such  investigations 
would  be  instructive  both  with  regard  to  predictions  from  theory  and  to  a judgment  of  the 
performance  of  the  shrouded  unit  as  a noise  source  in  comparison  with  an  unshrouded  pro- 
peller. 
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